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We present a regularized and renormalized version of the one-loop nonlinear relaxation equations 
that determine the non-equilibrium time evolution of a classical (constant) field coupled to its 
quantum fluctuations. We obtain a computational method in which the evaluation of divergent 
fluctuation integrals and the evaluation of the exact finite parts are cleanly separated so as to 
allow for a wide freedom in the choice of regularization and renormalization schemes. We use 
dimensional regularization here. Within the same formalism we analyze also the regularization 
and renormalization of the energy-momentum tensor. The energy density serves to monitor the 
' reliability of our numerical computation. The method is applied to the simple case of a scalar (f> 4 

theory; the results are similar to the ones found previously by other groups. 
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I. INTRODUCTION 



Phase transitions in elementary particle physics have been studied over the last two decades in various contexts. The 
phase transitions of grand unified theories may lead to inflationary periods of the early universe [Q; the electroweak 
phase transition can be based on a theory whose parameters - up to the Higgs sector - are meanwhile well known and 
studies both in perturbation theory and on the lattice show considerable progress in its understanding. The 
hadronic phase transition is being investigated both theoretically and experimentally || . 

For some applications it is sufficient to work in finite temperature quantum field theory, i. e. to consider only states 
\Q ' of thermal equilibrium. For the inflationary phase transition, at least, a real time study involving nonequilibrium 
ON ■ dynamics is necessary if one intends to describe the process of reheating. The simplest version of the inflationary 
scenario which consists in introducing a friction term [J7|,|| has been shown recently [|| to be on weak theoretical 
grounds. Linear and nonlinear relaxation have therefore been studied by various authors |p|-pT[. The one-loop 
relaxation equations studied here do not seem to lead to a proper thermalization, presumably they have to be modified 
in such a way as to include the interactions among the fluctuating fields. However, they will certainly describe an 
< initial stage of a more involved dynamical development. They present therefore a first step that has to be studied 
and well understood. Another subject to which nonequilibrium dynamics has been applied recently is the hadronic 
phase transition jl2| and the possible formation of disordered chiral condensates (l^ jl4|] . 

The computations of the nonlinear relaxation equations requires regularization and renormalization. Thus far 
P,|T4| noncovariant momentum cutoffs were used in such computations. While it would certainly not be an essential 
difficulty to replace those by a Pauli-Villars type regularization one may wish, in computations involving nonabelian 
gauge theories, to be able to use dimensional regularization as well. Furthermore it is more convenient to compute 
only convergent integrals instead of dealing with divergent integrals by varying their cutoff. 

It is the aim of this work to present a computational method in which the process of regularization and renormal- 
ization is cleanly separated from the numerical computation of finite parts of the one-loop integrals. The method is 
similar to the one developed in |15| which has been applied for the computation of one loop contributions to reaction 
rates associated with sphalerons |16[ , instantons [|l7| and bounces [|l8| . While in the cases mentioned the computations 
involved Euclidean Green functions which are behaved smoothly, here we will have to compute the trace of a Green 
function in Minkowski space, i. e. involving real time. The oscillating behaviour of such Green functions could present 
a new difficulty for the application of the computational scheme. 
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In this first presentation of the method we have studied its practical numerical application only in the simplest 
possible model, a scalar field theory with a A<£> self interaction. This theory has also been the subject of previous 
studies U and the results are similar. Especially it is found that the parametric resonance implied by an oscillating 
background field does not fully develop; back reaction damps the classical amplitude and the system relaxes to a new, 
essentially stationary oscillating phase. 

The plan of this work is as follows: in section || we recall the basic definitions and relations; in section III 



we 



present the one-loop nonlinear relaxation equations; we prepare the regularization in section IV by expanding the 
fluctuation modes in orders of the vertex function governed by the classical field and by deriving the large momentum 
behaviour of the first terms; regularization is then straightforward, the renormalization requires some algebra, both 
are presented in section ^; the formalism is extended , in section [vj| to the renormalization of the energy-momentum 
tensor; the numerical computation is discussed in section IVIJ; we conclude in section VIIl with a discussion of the 



numerical results and an outlook to more realistic and more general applications of the method. 



II. BASIC RELATIONS 

We restrict our study to self-interacting scalar 4 -theory without spontaneous symmetry breaking. The Lagrangean 
density is given by 

c = \d^d^ - i™ 2 $ 2 - ^<& 4 . (i) 

We split the field $ into its expectation value <p and the quantum fluctuations ip: 

®{x,t) = 4>{t)+^{x,t) (2) 

with 

0(t) = i)> = ^|) (3) 

where p(t) is the density matrix of the system which satifies the Liouville equation 

.dp(t) 



dt 



= [H(t),p(t)]. (4) 



The Lagrangean then takes the form 



C = Co + Cj (5) 



with 



1 

!„,„.., 1 o ,o A 



|5 M ^--m 2 2 --^ 4 , (6) 
A = drfd»><j> - m 2 V0 - ^ 4 - ^V 3 - \^ 2 4> 2 ~ ■ (7) 

This decomposition will be used as the basis for a perturbative expansion. Though we will consider the so-called one- 
loop equations which are nonperturbative, we will need such a perturbative scheme in order to define the one-loop 
summation and its renormalization. Since we are dealing with a nonequilibrium system the appropriate expansion is 
the CTP formalism (see e.g. p9|). We will consider the zero temperature case only. 
The Green function is defined by a 2 x 2 matrix 

G++(t,x:t',x') G + -(t,x;t',x') 

where 
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iG ++ (t,x;t',x") = (T<S>(t,x)<£(t,x)) (9) 

-iG—(t,2;t\tf ) = (T$(t,g)®(t',£)) (10) 

iG + -(t,x;t',x' ) = (11) 

iG~ + (M; £',£') = iG + -{t',x";t,x) . (12) 



Here the brackets (} refer to an expectation value with respect to the density matrix p. 
The Green functions can be decomposed as 

G ++ (t, x; t', x") = G>(t, x- 1', x')Q(t - t') + G< (t, x; t' , x')Q(t' - t) (13) 

G (t, x; t', x") = G>(t, x; t' , x')Q{t' -t) + G<(t, x; t' , x')Q(t - t') (14) 

G + -(t,x;t',x') = -G < (t,x;t',x") (15) 

G- + {t,x;t',tf) = -G > (t,x;t' 1 x') = -G < {t' ,x';t,x) (16) 

with 

G > (t,x;t',x l ) = i($(t,i)$(t',f )) . (17) 
These equations apply to the exact Green functions. For the free Green functions G > (t, x; t', a?) is given explicitly 

by 

/d 3 k i 

p-exp(-H(f-f)+ifc(x^1) (18) 



with ujk = \/k 2 + to 2 . 

The vertices are obtained from C\ ordering the terms in powers of ip 

Ci = {d^tj, - m 2 <jy - ^ 3 )V - " ^ " ^ 4 ■ (19) 

Every term corresponds to a vertex. In the matrix form of the CTP formalism the vertex operators are given by 

_f-*r n (t,2) o 



The subscript n denotes the number of fluctuation fields ip entering the vertex. We have 



tT 1 =Hp-m a -^{t))4>(t) (21) 

*r 2 = -i±<p (22) 

iT 3 = -i\(j> (23) 

iTi = -iX . (24) 



III. EQUATIONS OF MOTION 

The equation of motion for the field <f>(t) follows from the tadpole condition j20| 

(iP ± (x,t))=Q. (25) 
After a straightforward calculation one obtains Q the equation 

4>(t) + m 2 0(t) + ^{t) + I^( t )G++(0) = (26) 

which is represented graphically in Fig. [I]. Here G ++ is the H — h matrix element of the exact Green function in the 
background field <f)(t). It can be expanded perturbatively as 
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G — Go — Go^Gq + Gq^Go^Go 



G (l + r 2 G ) 



G ++ (t,x) satisfies the differential equation: 



dt 2 ~ 2 
So we have effectively a time-dependent mass term 



- A + m 2 + ^ 2 (i) ) G++(x; x') = (x - x>) 



(27) 



m 2 (t) = m 2 + ~<j> 2 {t) . 
Using translation invariance in x we introduce the Fourier transform 



(28) 



G£ + (t,t')= / d A x G ++ (t,x;t',0)eM-ikx) 



and denote the associated energy as 



We make the ansatz for G^ + (t, t') 



ui k (t) 



k 2 + m 2 + ^ 2 (t) 



Gt + (t,t') = C- 1 {U+{t)U^(t')6(t - f ) + U+(1?)Uj-(t)dtf t)} 
where the Uj^ are solutions of the homogenous problem and C _1 is the Wronskian. The initial conditions are 

^ + (0) = 1 C4 + (0) = -iu° k 
[/ fc "(0) = l l/fc(0)=to;g 



with 



fc 2 +m 2 + ^ 2 (0) 



(29) 



(30) 



(31) 



(32) 



(33) 



The complete Green function at equal times reads then 



Gt + (t,t) = ^\U+(t)\ 2 , 

where we have used U^(t) = U£(t)* . 

The resulting equation of motion for the classical field <f>(t) and the mode functions U£(t) are then 



(34) 



4 > (t)+m 2 ^t) + ^ 3 (t) + ^ ( j ) {t) 



d 3 k \U+(t)\ 2 



(2tt)3 2u,° 



We denote the fluctuation integral in (|35[) as 



U+(0) = 1 U+(0) = -iu° k 



A f d 3 fc \U+(t)f 



. . ,o i v A /■ d J , 



• 3 2 W ° 



(35) 

(36) 
(37) 

(38) 



It determines the back- reaction of the fluctuations onto the classical field (f>(t). 

An important check for the consistency of our numerical analysis will be the conservation of energy. The energy 
density is given by 
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(39) 



Calculating the trace over the Hamiltonian we obtain 
£ = ^(t) + imV(t) + ^ 4 (*) 



(40) 



Using the equations of motion it is easy to see that the time derivative of the energy density vanishes. The equations 
obtained so far are yet formal since they contain divergent quantities. We will present their renormalized form in 
section [v| 



IV. PERTURB ATIVE EXPANSION 



In order to prepare the renormalized version of the equations given in the previous section we introduce a suitable 



expansion of the mode functions U k {t). Adding the term ^(j) 2 (0)Ujt (t) on both sides of the mode function equation 



it takes the form 



with 



uUt) = -v{t)u+{t) 



(41) 



V(t) :=-(^(t)-f(0)) 



4. 



k 2 +m 2 + -</> 2 (0) 



1/2 



Including the initial conditions (32) the mode functions satisfy the equivalent integral equation 

oo 

U+(t) = e-<t + J dt'A ktIet (t - t')V(t')U+(t') 



with 



A fc , rc t(i - t') = *Q(t - t') sm(u° k (t - t')) 



(42) 



(43) 



(44) 



We separate U^(t) into the trivial part corresponding to the case V(t) — and a function fk(t) which represents the 
reaction to the potential by making the ansatz 



U+(t) = e-< t (l + Mt)) 



(45) 



fk(t) satisfies then the integral equation 



f k {t) = /dt'A fe , ret (i - t')V{t'){\ + / fe (0)e^ (t -° 



and an equivalent differential equation 

f k (t)-2iu°J k (t) = -V(t)(l + f k (t)) 
with the initial conditions /fc(0) = fk(Q) = 0. 
We expand now fk(t) with respect to orders in V(t) by writing 



(46) 



(47) 
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fk(t)=fi 1) (t) + f^{t)+fi 3 > (<) + .... (48) 
= /i 1) W + /f (*) (49) 



where f k (t) is of n'th order in V(i) and f^ n \t) is the sum over all orders beginning with the n'th one. The are 
obtained by iterating the integral equation ( j46| ) or the differential equation (f47|). The function f k (t) is identical to 
the function fk (t) itself which is obtained by solving (|47|) , the function f k (t) can be obtained as 

t 

ff\t) = / dt'A ktIet (t - OnO/FW^* -0 (50) 



or by solving the inhomogeneous differential equation 

/f>(t) - 2t W g/f>(t) = -V(i)/F(*) • (51) 

Note that in this way one avoids the computation of fj?\t) via the small difference fk(t) — This feature is 

especially important if deeper subtractions are required as in the case of fermion fields. 

The order on the potential V(t) will determine the behaviour of the functions fj^ at large momentum. We will 
give here the relevant leading terms for fjp(t) and f^ 2 \t). We have 

t 

/F(*) = £o / dt'(exp(2iu,°(i - if)) l)V(t') . (52) 



J k 



Integrating by parts we obtain 



t t 

i 



f * ){t) = -2J>J dt ' m ~ + iRF J ~ t')Mt') (53) 



or, by another integration by parts 



t 

i 



o 
t 



8( W ; 



0V3 



di' exp(2iu;£(t - t'))V(t') . (55) 



(2) 

Similarly we find for the leading behaviour of fj, (t) 



>k 

t t' 

/fc 2) (*) = -I7^m2 f dt ' fdt"V(t')V(t")+O((u; k )- s ). 



4(^) 2 



o o 



The leading terms of fjp(t) and /£ (i) in this expansion in powers of (w°) _1 are the same as for fjp(t) and fj^\t) 
respectively. 

Unlike a WKB ansatz for U^(t) the expansion presented here can be easily extendend Jll| to coupled channel 
systems and higher orders in the expansion (if deeper subtractions are required) . 



V. RENORMALIZATION 



The fluctuation term (38) occuring in the equation of motion (|35|) can then be written as 



G 



_ A f d 3 fc \e-^(l + f k (t))\< 



Inserting our expansion we obtain 

A f d 3 k 1 



(58) 



A f d 3 k 1 



2 7 (2tt) 



The second integral is convergent. Though Re/^ (i) and /^(t)! 2 behave separately only as (w°) 2 one finds, 
using Eqs. (p3|) and (p6|), that this leading behaviour cancels among the two contributions. This integral can be 



computed numerically using the functions and ff® obtained by solving (|4?| ) and (l5ll). In the first integral we find 
a quadratic divergence which corresponds to the tadpole graph in 4 -theory and a logarithmic one associated with 
2Rejf(i). 

Indeed the real part of fjp(t) is obtained from ([33]) as 

2Re/W(t) = --lj|nt)- /dt'co S (2^(t-t'))^(*')l • (59) 



o 



The first term behaves as (w°) -2 and leads therefore to a logarithmic divergence in AM 2 . The second one yields a 
finite contribution. It will be included into the finite part of AAA 2 . 
So altogether we have 

a * > 2 / \ A f d 3 k 1 A f d 3 k 1 V(i) . , .n , , 



2 7 (27T) 3 2^o 2 7 (2tt) 3 2^° 2(c°) 2 
with 

t 





-L{ 2 Re/f (i) + |/«(*)| 2 } • (61) 



A 

+ 2 J (2tt) 3 2u/ 

The first term in (|6^) corresponds to the zeroth order diagram in Fig. |l|. The second term in (|6C| ) and the first term 
in (|6l]) together correspond to the first order diagram in Fig. |J the second term in ( |6l| ) sums up all the higher order 
diagrams. By this term-to-term equivalence of our expressions to diagrams of CTP perturbation theory we have 
achieved here a clean separation between finite quantities that can be computed numerically and the renormalization 
parts that are computed analytically up to Fourier transforms of the external sources. We consider this to be an 
important feature of our method, as a the regularization may now be chosen freely, as required in order to maintain 
the symmetries of the theory. 



The first divergent integral can be rewritten as 

d 3 k 1 f d 4 fc 



(2tt) 3 2uj° k J ( 2 tt) 4 k 2 - to 2 + to ' 



(62) 



an identity that holds in dimensional as well as in Pauli-Villars regularization. The integral is the one associated to 
the tadpole graph Fig. |^a and can be absorbed into the renormalization of the mass term. In the same way the second 
integral is equivalent to 

d 3 fc 1 If d 4 fc -i 

(63) 



(2tt) 3 8(w°) 3 2 J (2tt) 4 (fc 2 - to 2 + io) 2 ■ 

It is associated with the Feynman graph Fig. |2|b which leads to a coupling constant renormalization. We will use 
dimensional regularization and renormalize the four point function at vanishing external momenta. We first rewrite 
the basic equation of motion, including appropriate counter terms, as 
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^ 2 1 ■ A + — ^ 3 (t) + A.M 2 (^(t)=0. 



0(t) + (m 2 + <Jm^)0(t) + 



(64) 



Next we separate from AAi 2 (t)(p(t) dimcnsionally regularized divergent terms 



d 3 ^ i 



(2tt)3 2u,o 



rcg 



\ml<f)(t) 
32tt 2 



d d fc 
(27r) d 



fc 2 + m 2 + ^02(0) 



In 



47T/T 
m 2 , 



-7 + 1 



(65) 



and 



A 



M)V(t) 



d 3 /c 



(2^) 3 4u;f 



(^) £ A0(t)V(t) 



A^(i)y(i) f 2 
32^2 \7 



d d /c 
(2^7 



fc 2 +m 2 + |^<?!> 2 (0) 



- + In 



47r^ 2 



(66) 



where we have introduced 



m 2 : =m 2 + - M £ </> 2 (0) . 
Recalling that V(t) = (\/2)(<p 2 (t) — 4> 2 {0j) these two expressions combine into 



(67) 



- 5m 2 (j>{t) 



(5A/z e 



A 



■ In ■ 



32tt 2 ml 



m 2 cb(t) + ^ 3 (t) 



64tt 2 



(68) 



where we have defined the rcnormalization constants 



r 2 Am 2 f 2 , 47T/Z 2 



7 + 1 



(69) 



and 



3A 2 (2 , 47ru 2 
- +ln — V- 



32tt 



m^ 



'7 



Note that these counterterms are independent of m and therefore from the initial value of the field </>(0) 
divergent integrals do depend on the initial condition. 
The renormalized equation of motion now reads 

4>{t) + (m 2 + Am^(t) + *±^V (t) + AMUtm) = 



(70) 
though the 

(71) 



with the finite corrections 



Am 2 - 



Am m 
o In — o 

32tt 2 m 2 



3A 2 , m 
AA = --r-^ln 



64tt 2 



2 (0) 



32tt 2 m 2 ' 

This equation is in a form well suited for numerical computation. 



(72) 
(73) 



VI. THE ENERGY-MOMENTUM TENSOR 



The renormalization scheme applied so far concerned the Green functions of the theory. The energy density on the 
other hand is part of the energy-momentum tensor, ft contains additional divergent terms which have to be defined by 
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a regularization and to be removed by new counter terms. These counter terms have been considered in the literature 
[|2 1| — |2-4J| . If one considers a space with vanishing curvature one has only two possible counter terms besides the ones 
already introduced into the Lagrangean so that 



T™ = T£? + 9lw ( SA + \8m 2 4> 2 {x) + ^(x) ) - .U ;ll! ,.0„<r - i),AA<ru) . (74) 



In the case considered here the second term does not contribute to the energy density. It will be useful, however, to 
consider the entire energy-momentum tensor including the space components 

Tij = -P9ij (75) 
where p is the pressure given formally (i. e. without regularization) by 

p = m + /H ^ (i W + 1 wci") - £ ■ <™> 



The constants A and SA have been determined by various authors |21|-|24|] using different regularizations. Coleman 
and Jackiw use the covariant Pauli-Villars regularization and Christensen the point-splitting technique and heat kernel 
regularization. Cooper et al. find a way of dealing with the problems of a three-momentum cutoff. We will use the 
covariant dimensional regularization as we did above in regularizing the fluctuation integral. 

In order to rcnormalize the energy we introduce the available counter terms into the unrenormalized expression 
( ffO| ) so that it reads now 



d 3 k 



^(l + 2Re/«W + |/i 1) (t)| 2 ) 



1 fW/ 



2 



(2tt) 3 [ 4 

+-M + + 2Re/F(t) + l/FWI 2 )} • (77) 

While SX and 5m 2 have already been fixed above we have to determine now the 'cosmological constant' SA; in order 
to do so we analyse the quartically divergent part of the fluctuation integral 

_ f d 3 fc 

tquart -J (2?r) 3 2 ■ U»J 

In dimensional regularization it becomes 

ml f 2 4tt^ 2 3 



^---64^|7 + ln ^ + 2--r 

This expression obviously depends on the initial condition via the 'initial mass' too. If we insert the expression (|67j ) 
we find 

^ uart = -^(m 4 + ATOV(0) + ^ 4 (0))^ + ln^ + ^- 7 ) . (80) 



The divergent parts proportional to 2 (O) and ^ 4 (0) will be cancelled by divergent parts of the fluctuation integral 
(see below). We define here the 'cosmological constant' counter term by 



in /2 4nu 2 3 \ 

6A =64^ 7 +ln ^- 7+ 2) ' (81) 



again independent of the initial condition, and retain in the expression for the energy a finite term 







647T 2 Ttl V ' 

Having considered the quartically divergent term in the fluctuation integral we turn now to the quadratic ones. 
One can show that the terms uj^Refj^\t)/2 and — Re(if^*(t))/2 cancel. The term proportional to gives in 

k 

dimensional regularization 



-'quad 



d 3 fc V(t) AmL,,,, „„, 5 , (2 , 4?ru 2 A 



(2?r) 3 AujI 64tt 



These terms have to be considered together with the logarithmically divergent ones which are proprotional to V(t) 2 . 
Their explicit form follows from an analysis of the leading behaviour. It combines to 

£ ^=- v ^Sw?^'-^ m2 - m2)2 £ + *%-^ ' (84) 

Instead of removing the complete terms \fP~' (t)| 2 /4w° and 2V(i)Re/ ) j; 1 ' ) (i)/4w° from the fluctuation integral we just 
subtract the leading behaviour from the integrand, this is algebraically less cumbersome at the expense of inducing 
a small difference of large quantities (which was found to be tolerable). As before (see (|68|)) the quadratically 
and logarithmically divergent terms combine in such a way that the mo is replaced by m in the divergent parts. The 
terms proportional to (j) 2 (t) and <j) A (t) are cancelled by the mass renormalzation and coupling constant renormalization 
counter terms up to finite contributions. The terms proportional to (f> 2 (0) and <f> (0) are cancelled by the corresponding 
terms occuring the quartically divergent part ([79]), again up to a finite remainder 

AA' = ^ (y^ 4 (0) - AmV(O)) . (85) 
We finally obtain for the energy density the renormalized expression 

Sren = \tf{t) + \{™ 2 + Am 2 )0 2 (t) + *±^*( t ) + AA + AA' 



^ ( 2Re/f(t) + |/«(*)| 2 ) 



d 3 fc 

(2tt) 3 [ 2 

+ 4^'^' 2 

+ ^(2Re/FW + l/FWI 2 ) + 0}- (86) 

In the expression ( f76|) for the pressure the energy £ has to be replaced by the renormalized one. This replacement 
absorbs all the counter terms proportional to (see ([74])). The remaining counter term has to be added and the 
constant A has to be chosen so as to cancel the remaining divergencies. Replacing at the same time the mode functions 
U£ by their expansion in terms of the we find 

d 2 

Prcn = -£rcn + 4> 2 (t) + A — (j) 2 {t) + pfl uct (87) 



dt 2 



with 

Pfluct 



+ |/F(<)| 3 - 2u> (l + /«(«)) /f>*(t)} . 



This expression looks at first sight hopelessly divergent and there is just one counter term to cancel these divergencies. 
Using dimensional regularization we find that the leading quartic divergence reduces to a finite term via 
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d 3 k 1 „ n ,2 k 2 . 



(2tt) 3 2ul 



The terms leading to a quadratic divergence are 



rog 



III , 



96tt 2 ' 



(89) 



u;°Re/«(t) Re(if^*(t)) + -^Ref^(t) 

oLOu 



(90) 



Analyzing their leading behaviour one finds that the possible quadratic divergence reduces again to a finite result via 



f d 3 k 1 f d 3 k 1 



Win 



V(t)ml 
48tt 2 



After evaluating these leading divergencies we find 
m\ V{t)ml 



7 + 1 )+V(t) 



48tt 2 



- + In 



ran 



(91) 



Pfluct 



96tt 2 48tt 2 



+ 



d 3 k 1 



((^) 2 + 



4 K) 7 



di sin2w^(t-i / )V r (i') 



+2Re(-^/( 1 )(t)-^/^^(t)/^*(i)) 







(92) 



The analysis of the logarithmically divergent terms becomes now very cumbersome. After some algebra and integra- 
tions by parts one finds the logarithmically divergent term 



, , d 3 fc 1 



(2tt) 3 2ujI 1 8(w 0) 



1 W + £ 



1 



4Kr 



rog 



v(t) 

96tt 2 



, , 2 4 V 
/U - + In =- - 7 



(93) 



This determines the divergent part of the renormalization constant A. The finite part can be chosen freely and 
constitutes a free parameter of the theory, used e.g. for the 'improved' energy-momentum tensor. We dispose of this 
freedom by choosing 



A = 



X 



1927T 2 



2 , 4tt/x 2 



(94) 



since then the expression in parentheses is just the standard equivalent of ln(A 2 /m 2 ) in Pauli-Villars regularization. 
The rcnormalizcd pressure then takes the final form 



Pron = -£ren + 9 (t) 



< v(t) 

96tt 2 48tt 



2"'0 



V(t) 
288tt 2 



+ 



r d 3 fc i 
J 



(^) 2 + ^)(2Re / f W + |/«(t)| 2 ) 
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t 



' "'" 1 1 dt cos2u)° k (t - t') V (t') 



6(ou°) 2 24( W «) 4 
1 ml 



cos (2u° k t)V(0) 



Vl2(^) 2 24(^) 4 , 
H/F(*)r-2^Re(»/?>(t)+i/«(t)/W*(t))L (95) 



VII. NUMERICAL ANALYSIS 

As discussed in section 3 the computation of the fluctuation term AM. 2 (t) which determines the back- reaction 
of the quantum fluctuations to the classical field <p(t) can be reduced to the computation of the mode functions 



/£ ' (t) — fk(t) and f k (t) via (|6l|). The differential equations satisfied by these functions, ([471) and (51), have been 



solved using an improved Runge-Kutta-scheme with six steps |25|] . We have also checked the accuracy of (t) by 
computing it via the alternative expression (|5C|). 

The momentum integrals are finite; however, the convergence is not very good, as the integrand decreases only 
as an inverse power of u> k . Furthermore the integrand is oscillating rapidly. If we denote the upper limit of the 
momentum integration as K the integrals occuring in AM 2 , the energy £ and the pressure p converge only as 
I(K) = loo + C ' /{<jj° K ) 2 to their value /qq. We have determined the values 1^ by fitting the integrals as a function of 
K to a function of this type. This way the numerical integration can be limited to moderate values of K ~ 20 -f- 40. 
A delicate part of the momentum integration is the region of small momenta. There the mode functions develop a 
parametric resonance before their back-reaction to the classical field (f>(t) suppresses the amplitude of the latter one 
and the system "shuts off" jlO]], i. e. reaches a stationary oscillating behaviour. We have chosen momentum intervals 
varying from Ak ~ 10 -4 at small momenta to Ak ~ .1 for larger ones. The time steps were chosen as At ~ 0.001. 
The quality of the numerical procedure was monitored by the conservation of energy which was fulfilled with an 
accuracy of better than 3%. With the accuracy given above we had to compute the mode functions for 1000 values of 
momentum. The computations can be performed on a modern PC or small workstation, they take then several hours 
of CPU time. 



VIII. RESULTS AND CONCLUSIONS 



We have presented here a computational scheme for solving the one-loop nonlinear relaxation equations with a 
manifestly covariant regularization and renormalization. 

The results are presented in Figs. 3 - 9 for two sets of parameters, similar to those of [||: X/8tt 2 — 0.1, m = 1 (as 
the general unit) and </>(0) = 5, <h(0) = and (f)(0) = 1, 0(0) = 0, respectively. In Figs. || and ||we display the classical 
amplitude (f>(t), in Figs. Q and [|we show the growth of the fluctuation integral AA4 2 (t) (see (|6l])). Fig. [s] shows the 
absolute value of the integrand (to be integrated with the measure k 2 dk/u> k ) of the second fluctuation integral in ( |6l| ) 
vs. u>u in a double-logarithmic scale for the same set of parameters and for t = 37.5. The amplitude of the integrand 
is seen to fall off as (u) k )~ A . In Fig. || we plot the classical and fluctuation energies as well as the total energy for 
the first parameter set. The pressure is displayed for the same parameter set in Fig. 0. Its average as t — > oo is seen 
to be £ /3 as found previously in ]23j . These numerical results show the same qualitative and similar quantitative 
features as those obtained in ^] . An initial large amplitude oscillation leads after a short time interval - depending on 
the initial amplitude - to a strong excitation of the low momentum fluctuation modes as expected from a parametric 
resonance. Then after a considerable decrease of the classical amplitude an essentially stationary oscillating regime is 
reached. In contrast to a computation in Hartree-approximation in Q we find that even for (f>(0) = 1 a considerable 
amount of energy is transferred to the fluctuation modes. The main difference between a small and large initial value 
of the classical field consists in the time that is needed for the fluctuations to build up. 

The field theory we have considered here was a rather simple one; as discussed in P,[To| one would not expect the 
one-loop approximation to be a good here nor has the system been found previously to display thermalization. Here 
this simple model served as a toy example for demonstrating the method. More interesting systems are of course 
gauge theories with many light degrees of freedom. There the differential equations for the gauge field fluctuations 
form coupled systems. We should like to stress that the method as developed in sections IV and ^ can be easily 
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generalized to deal with such coupled channel systems. This generalization is entirely analoguous to the treatment of 
coupled channel one loop computations in 
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FIGURE CAPTIONS 



FIG. 1. Diagrammatic representation of the one- loop equation (|26|). 

FIG 2. Renormalization parts: a) Tadpole diagram, Eq. (|62[), b) Fish diagram, Eq. (^). 

FIG 3. The classical amplitude <j>(t) for A/8tt 2 = 0.1 and 0(0) = 5. 

FIG 4. The fluctuation integral AM 2 (t) for the same set of parameters. 

FIG 5. The fluctuation integrand as explained in the text. The absolute value of the integrand is plotted vs. on a double 
logarithmic scale. The straight line indicates a power behaviour as (w')" 4 . 

FIG 6. The classical and fluctuation energies as a function of time for the same parameters: classical energy (short dashed 
line), the fluctuation energy (long dashed line) and total energy (solid line). 



FIG 7. The pressure as a function of time for the same parameter set. The horizontal line indicates the value p — £/3. 
FIG 8. The classical amplitude 0(i) for A/8tt 2 = 0.1 and 0(0) = 1. 
FIG 9. The fluctuation integral AM 2 (t) for the same set of parameters. 
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